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THE STELLAR VELOCITY DISPERSION OF THE LENS GALAXY IN MG2016+112 AT Z=1.004i 
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ABSTRACT 

We present a direct measurement of the stellar velocity dispersion of the early- type lens galaxy D in the 
system MG2016-I-112 {z = 1.004), determined from a spectrum obtained with the Echelle Spectrograph 
and Imager (ESI) on the W.M. Keck-II Telescope, as part of the Lenses Structure and Dynamics (LSD) 
Survey. We find a velocity dispersion of a^p — 304 ± 27 kms^^ inside an effective circular aperture with 
a radius of 0.65", corresponding to a central velocity dispersion of cr = 328 ± 32 kms~^. From a Hubble 
Space Telescope F160W-band image, we measure the effective radius and effective surface brightness 
in order to determine the offset of the lens galaxy with respect to the local Fundamental Plane. The 
offset corresponds to an evolution of the rest-frame effective mass-to-light ratio of A\og{M/ Lb) = 
—0.62 ± 0.08 from z = to z = 1.004. By interpreting colors and offset of the FP with two independent 
stellar population synthesis models, we obtain a single-burst equivalent age of 2.8±0.8 Gyr (i.e. Zf>1.9) 
and a supersolar metallicity of log[^/.Z0]=O.25±O.25. The lens galaxy is therefore a massive elliptical 
dominated by an old and metal rich stellar population at z > 1. The excellent agreement of the stellar 
velocity dispersion with that predicted from recent lens models confirms that the angular separation of 
the multiple images of the background QSO is predominantly due to the lens galaxy, and not to a massive 
"dark cluster" , in agreement with recent weak lensing and X-ray observations. However, the significant 
overdensity of galaxies in the field might indicate that this system is a proto-cluster, in formation around 
galaxy D, responsible for the '--^10% external shear inferred from the strong lens models. 

Subject headings: gravitational lensing — galaxies: elliptical and lenticular, cD — galaxies: evolution 
— - galaxies: formation — galaxies: structure 



1. INTRODUCTION 

In the cold dark matter (CDM) cosmological scenario, 
structures in the Universe form through hierarchical merg- 
ing of smaller structures (White & Rees 1978; Blumenthal 
et al. 1984; Davis et al. 1985). Within this general frame- 
work, early-type galaxies (E/SO) in the cores of rich clus- 
ters form at high redshifts (z > 2), corresponding with the 
first dark-matter overdensity peaks, as opposed to a later 
formation epoch for field E/SO (Kauffmann 1996). Clus- 
ters subsequently form around these seeds by accretion 
of smaller-mass galaxies, with significant structural and 
dynamical evolution occurring at more recent cosmologi- 
cal times. Recent observations have shown that massive 
cluster E/SO were already assembled at z > 1 and sub- 
sequently evolved passively through ageing of their stellar 
populations (e.g. van Dokkum et al. 1998; Stanford, Eisen- 
hardt & Dickinson 1998). Similarly, field E/SO galaxies 
seem not to evolve dramatically between z = 1 and z = 0, 
both in number (Schade et al. 1999, Im et al. 2002) and 
structural properties (van Dokkum et al. 2001; Treu et al. 
1999, 2001b, hereafter TOlb), although secondary episodes 
of star formation might be frequent at z < 1 (Menanteau 
et al. 2001; Treu et al. 2002, hereafter T02). However, 
most observational results on the evolution of E/SO con- 
cern the evolution of their stellar populations and little is 



known about the evolution of their internal structure. 

To comprehensively quantify the luminosity, color and 
structural evolution of the stellar mass and of the dark- 
matter halo of E/SO as function of redshift, we are con- 
ducting an observational program with the Echelle Spec- 
trograph and Imager (ESI) on the W.M. Keck-II Tele- 
scope: the Lenses Structure and Dynamics (LSD) Survey. 
Aim of the LSD survey is to measure the internal kine- 
matics of a dozen gravitational- lens galaxies up to z = 1, 
allowing a powerful combination of dynamical and lensing 
constraints on their mass distribution. The LSD Survey 
its main goals and observing strategy will be described in 
detail elsewhere (Treu & Koopmans 2002, in preparation). 

Here we present the first result of the LSD Survey, a 
measurement of the stellar velocity dispersion of the lens 
galaxy in the system MG2016-I-112 at z = 1.004. A sum- 
mary of relevant prior observations and a new lens model 
can be found in Koopmans et al. (2002) and references 
therein. The primary lens galaxy (D) in MG2016-I-112 
is the highest spectroscopically-confirmed redshift lens 
galaxy known to date. The suggestion that this lens was 
embedded in a massive "dark cluster", based on ASCA 
X-Ray observations (Hattori et al. 1997), was recently 
shown to be incorrect by high-resolution Chandra obser- 
vations, showing no evidence for hot X-ray gas (Chartas 



Based on observations collected at W. M. Keck Observatory, which is operated jointly by the California Institute of Technology and the 
University of California, and with the NASA/ESA Hubble Space Telescope, obtained at STScI, which is operated by AURA, under NASA 
contract NAS5-26555. 

^ see also http://www.its.caltech.edu/~koopmans/LSD or http://www.astro.caltech.edu/~tt/LSD 
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et al. 2001). On the other hand, deep optical studies 
that show an overdensity of at least six bright E/SO with 
the same redshift as the lens galaxy (Benitez et al. 1999; 
Soucail et al. 2001; Clowe et al. 2001). The absence of 
a significant weak-lensing signal (Clowe et al. 2001) and 
the absence of X-ray emission, however, shows that these 
galaxies are not associated with a massive evolved cluster. 

In the following, we assume for definiteness that the 
Hubble constant, the matter density, and the cosmologi- 
cal constant are Hq = 65 kms^^Mpc^^, Jim — 0.3, and 
fl\ = 0.7, respectively. 



2. OBSERVATIONS 

2.1. Keck Spectroscopy 

We observed MG2016+112 using ESI on the 
W.M. Keck-II Telescope during four consecutive nights 
(23-26 July, 2001), with a total integration time of 8.5 hrs. 
The seeing was good (0'.'6-0'.'8) and three out of four nights 
were photometric. Between each exposure of 1800 s, we 
dithered along the slit to allow for a better removal of sky 
residuals in the red end of the spectrum. The slit (20" in 
length) was positioned at a position angle of 13° in order 
to include galaxies D and E, and arc C (see Fig. 2). 

The slit width of 1'.'25 yields an instrumental resolu- 
tion of 30 km which is adequate for measuring the 
stellar velocity dispersion and remove narrow sky emis- 
sion lines. Data reduction was performed using a newly 
developed IRAF package^, that combines the procedures 
described in Treu et al. (1999, 2001a) and Smith et al. 
(2001) for the treatment of echelle distortions and sky 
subtraction. Spectra of G-K giants observed at twilight 
with a 0.3" slit width were used as stellar templates, after 
appropriate smoothing to match the instrumental resolu- 
tion of the 1.25" slit. The best fit was obtained for a 
G4-III stellar template, using the Gauss-Hermite pixel- 
fitting software on the spectral region around the G band 
(Fig. 1). This yielded a stellar velocity dispersion of 
Cap = 304 ±27 km s~^ inside an effective circular aperture 
with radius 0"65. The scatter in the fits for the different 
stellar templates was used to estimate the uncertainties re- 
lated to template mismatches (see TOla). The error on Cap 
therefore includes both the random error contributions and 
the systematic uncertainties due to template mismatches, 
continuum fitting, and fitted spectral range. The inferred 
central velocity dispersion is tr = (1.08 ± 0.05) • Cap, i. e. 
fj = 328 ± 32 kms~^. The correction is based on the aver- 
age ratio (and spread) between the central stellar velocity 
dispersion (i.e. measured in a Re/S aperture) and the dis- 
persion measured from a larger aperture (~ 2i?e in our 
case), as determined from local E/SO (see TOla for de- 
tails). The predicted velocity dispersion from isothermal 
gravitational- lens models, i.e. 320-340 km s^^ (Koopmans 
et al. 2002), is in excellent agreement with the spectro- 
scopic determination of the central stellar velocity disper- 
sion (see also Kochanek 1993, 1994 and Kochanek et al. 
2000 for a general discussion). In addition, we detect [Oil] 
emission from galaxies D and E at nearly identical red- 
shifts (z = 1.004; Fig.2). 

2 EASI2D, developed by D. Sand and T. Treu; Sand et al. (2002), in 
^ Obtained as part of the CASTLeS Survey 
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Fig. 1. — The G-band absorption feature at ~4304 A in the spec- 
trum of galaxy D in MG2016-I-112. The rest-frame wavelength scale 
is shown on the top axis for reference. The smooth solid line is the 
best-fitting template (spectral-type G4— III) convolved to the mea- 
sured velocity dispersion (Tap = 304±27 kms~^. The dashed region 
is affected by sky-subtraction residuals and was omitted in the fit. 

2.2. Hubble Space Telescope Imaging 

Optical and infrared Hubble Space Telescope (HST) im- 
ages of the system are available from the HST archive^. 
The Near Infrared Camera and Multi Object Spectrograph 
(NICMOS) Camera 2 imaged the galaxy in F160W for 
5112s, while the Wide Field and Planetary Camera 2 im- 
aged the system through filters F555W and F814W for 
5200s each (Fig 2). The images were reduced and surface 
photometry performed on the F160W and F814W images 
as described in Treu et al. (1999, 2001a), except that cos- 
mic rays rejection on the F814W image was performed 
using the LA Cosmic Algorithm (van Dokkum 2001). The 
lens galaxy is very faint at F555W and barely detected in 
the Planetary Camera images, yielding a very uncertain 
flux. The galaxy brightness profiles are well represented 
by an R'^/* profile, but the F160W- and F814W-band ef- 
fective radii are significantly different and the galaxy shows 
a significant reddening towards its core. We attribute the 
larger effective radius in F814W-band to a bluer extended 
envelope around galaxy D (seen as a "wing" in the bright- 
ness profile), which presumably comes from a younger stel- 
lar population that has been accreted from smaller galaxies 
(e.g. galaxy E). Because the F160W-band hght is a better 
tracer of the underlying old stellar population (dominant 
in mass) and is not contaminated by the [Oil] emission 
line (Fig.2), we only use Re.Fieow for our analysis. In 
addition, the fits provided the integrated magnitudes and 
effective surface brightnesses. The relevant observational 
quantities of MG2016-I-112 are listed in Tab. 1. Rest frame 
photometric quantities are corrected for galactic extinction 
using E(B-V)=0.235 from Schlegel et al. (1998). 
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Fig. 2.— HST images of MG2016+112: WFPC2 F555W (V) and 
F814W (I) images, the NICMOS-2 F160W (H) image. The exposure 
times are 5200s (V and I), and 5120s (H), and the pixels scale are 
0.0455" (V), 0.1" (I), 0.075" (H). The two lensed quasar images (A 
and B), the primary lens galaxy (D), the arc C and the additional 
galaxy (E) are indicated on the H-band image. The emission from 
the [Oil] doublet detected from galaxies E and D is shown in the 
bottom right panel. Note that the doublet in D is blended by the 
large velocity dispersion. 

3. RESULTS 

Using the stellar velocity dispersion, effective radius, ef- 
fective surface brightness and colors of the lens galaxy, wc 
can now quantify the evolution of its stellar mass-to-light 
ratio, as well as the age and metallicity of the dominant 
stellar population. A discussion is given in Sec. 4. 

Recent studies have shown that E/SO in clusters and 
in the field define a tight Fundamental Plane (FP) out to 
zr^O.7- 0.8 (van Dokkum et al. 1998; T02), with slopes 
very similar to the ones observed in the local Universe. As- 
suming that galaxy D lies on a FP with slopes similar to 
those in the local Universe'*, we can obtain the evolution 
of the intercept 7 of the FP with redshift, which is re- 
lated to the evolution of the average effective mass to light 
ratio (M/L). In Fig. 3 we plot the evolution of the M/L 
for cluster and field E/SO as function of redshift (small 
squares and pentagons; from T02), together with the value 
obtained for galaxy D (large filled square). From the lo- 
cal FP, one finds an offset ASBe = —1.55 mag arcscc^, 
which corresponds to A log M /Lb = —0.62 ± 0.08 for pure 
luminosity evolution, the error being dominated by the ob- 
served FP parameters of galaxy D. We also note that errors 
on Re and SBg are correlated and partly cancel in the FP 
relation (e.g. Treu et al. 2001a; herafter TOla). The M/L 
evolution is intermediate between that of clusters and field 
E/SO, when we assume that a linear fit, Alog M/ Lb cx z 
(see Fig. 3), provides a good description of the data and 
can be extrapolated to z = 1.004. Eventhough this is 
only a single-galaxy measurement - not a determination 

* In the following, we adopt log iie + log/iso = 1-25 log cr + 0.32 SBg 
TOlb for discussion). The uncertainties on the local determination of 




of the FP at z ~ 1 - the inferred M/L evolution suggests 
that passive evolution of and old stellar population is a 
good representation of the evolution of massive E/SO out 
to 2; = 1.004. Additional measurements, covering a range 
of FP parameters, are needed to confirm this result. 
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Fig. 3. — Evolution of the M/L as inferred from the evolution of 
the FP. Filled pentagons and open squares represent field and clus- 
ter measurements, respectively, while the thick full and dashed lines 
indicate linear fits to the M/L evolution for field and cluster E/SO 
(sec T02 for details). The large filled square at 2: = 1.004 indicates 
the M/L evolution of galaxy D in MG2016-I-112. 



The color and M/L evolution of lens galaxy D can be 
used to infer age and metallicity of its stellar population. 
A single-burst stellar population, generated with two in- 
dependent population synthesis codes (Bruzual & Char- 
lot 1993; Fioc & Rocca-Volmerange 1997), gives a single- 
burst equivalent age of 2.8±0.8 Gyr (i.e. Zf >1.9) and a 
supersolar metallicity log[Z/Z0]=O.25±O.25 (Fig. 4), as- 
suming there is no significant internal extinction in the 
lens galaxy. The contours in Fig. 4 are based on the 
difference between the model and the observed color and 
mass-to- light ratio evolution (and their errors). Because 
the colors and the luminosity can be dominated by a small 
mass-fraction of bright young stars, this age should be re- 
garded as a lower limit for age of the dominant stellar 
population. The supersolar metallicity is slightly higher, 
but within 1-a consistent with results obtained from E/SO 
in clusters at redshifts z w 0.3 — 0.9, which have been as- 
sembled through hierarchical merging of smaller galaxies 
(e.g. Ferreras, Chariot & Silk 1999). 



4. DISCUSSION & CONCLUSIONS 

We have presented a direct measurement of the 
stellar velocity dispersion of the lens galaxy in the 
highest spectroscopically-confirmed redshift lens system 
MG2016-hll2 {z = 1.004), as part of the Lenses Struc- 
ture and Dynamics (LSD) Survey. The effective radius 
and surface brightness are determined from HST images, 
to compare the properties of the lens galaxy with those of 
E/SO in the local Universe. The colors and evolution of 
M/Lb of its stellar population, as derived by comparison 
to the local FP, indicate an old (~3 Gyr) metal-rich stel- 
lar population, consistent with a passively evolved single 

— 8.895 (Bender et al. 1998) for the local FP in rest-frame B-band (see 
the FP are negligible in our analysis. 



4 



Koopmans & Treu 



stellar population formed at Zf > 1.9. The large stellar ve- 
locity dispersion and metallicity of the lens galaxy is com- 
parable to the most massive present-day cluster galaxies 
(Ferreras et al. 1999; Fisher et al. 1995). In addition, 
the F814W-band surface brightness profile at large radii 
deviates slightly from the i?^'^ profile, possibly indicating 
an outer envelope, as typically found for brightest cluster 
galaxies and cD galaxies (Graham et al. 1996). This bluer 
envelope and the [Oil] emission that we detect from galaxy 
D (see also Soucail et al. 2001) could be due to the accre- 
tion of a younger stellar populations from smaller gas-rich 
galaxies, such as the closest companion to galaxy D for 
which we detect [Oil] in emission (Fig. 2). 

I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' 
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Fig. 4. — Likelihood contours of the age and metaUcity of the 
stellar population of MG2016+112, galaxy D. The inner and outer 
contours indicate the 68% and 95% probability. Panels (a) and (b) 
are for population synthesis models from Bruzual & Chariot (1993) 
and Fioc & Rocca-Volmerange (1997), respectively. 

This scenario is consistent with the gravitational-lens 
system being embedded in a galaxy overdensity (i.e. a 
proto-cluster; see also Benitez et al. 1999), which is not 
yet relaxed to a more centrally concentrated system and 
still shows evidence for ongoing merging and accretion. 
In addition, it would explain the absence of a significant 
weak-lensing signal (Clowe et al. 2001), the relatively small 
lensed-image separation, as well as constraints from re- 
cent Chandra X-ray observations (Chartas et al. 2001). 
In particular, the spread of the velocities of the galaxies 
in the overdensity (Soucail et al. 2001), if they are as- 



sumed to be virialized, implies a mass three times larger 
than the S-cr upper limit inferred from the non-detection 
of extended X-ray emission with Chandra (Chartas et al. 
2001). We therefore think that the assumption of virial- 
ization is not correct and that these galaxies form a (non- 
virialized) proto-cluster. Three galaxies in the field are 
spatially clustered dX z — 0.97 and might be falling to- 
wards the overdensity with high velocity (Soucail et al. 
2001), consistent with a cluster in formation. 



Redshift (D) 


1.004±0.001 


FimW (mag) 


18.24±0.02 


F&UW-FimW (mag) 


3.3±0.1 


SBe,_Fi60H' (mag/arcsec^) 


17.64±0.40 


Re,FsiiW (arcsec) 


0.65±0.10 


Re,Fiww (arcsec) 


0.31±0.06 


cr*(< 0.65") (kms-i) 


304±27 


b/a={l - e) 


0.75±0.10 


Major axis P. A. (°) 


121 ±2 


(7 (kms^^) 


328±32 


B — I (mag) 


1.98±0.15 


Mb (mag) 


-22.53±0.10 


Mi (mag) 


-24.51±0.04 


SBe^B (mag/arcsec^) 


18.12±0.50 


SBej (mag/arcsec^) 


16.14±0.40 



Table 1. — Observed spectro-photometric quantities of 
galaxy D in MG2016-)-112. The second part of the ta- 
ble lists rest- frame quantities, derived from the observed 
quantities (see text). 
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